Nonlinear Dynamics in the Resonance Lineshape of NbN Superconducting Resonators 
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In this work we report on unusual nonlinear dynamics measured in the resonance response of 
NbN superconducting microwave resonators. The nonlinear dynamics, occurring at relatively low 
input powers (2-4 orders of magnitude lower than Nb), and which include among others, jumps 
in the resonance lineshape, hysteresis loops changing direction and resonance frequency shift, are 
measured herein using varying input power, applied magnetic field, white noise and rapid frequency 
sweeps. Based on these measurement results, we consider a hypothesis according to which local 
heating of weak links forming at the boundaries of the NbN grains are responsible for the observed 
behavior, and we show that most of the experimental results are qualitatively consistent with such 
hypothesis. 

PACS numbers: 85.25.-j, 74.50.+r, 47.20.Ky, 05.45.-a 



I. INTRODUCTION 

Understanding the underlying mechanisms that cause 
and manifest nonlinear effects in superconductors has a 
significant implications for both basic science and tech- 
nology. Nonlinear effects in superconductors may be ex- 
ploited to demonstrate some important quantum phe- 
nomena in the microwave regime as was shown in Refs. 
PI Q, and as was suggested recently in Refs. 0, E|- 
Whereas technologically, these effects can play a positive 
or negative role depending on the application. On the one 
hand, they are very useful in a wide range of nonlinear 
devices such as amplifiers 0, @L mixers , single pho- 
ton detectors 0, and SQUIDs On the other hand, in 
other applications mainly in the telecommunication area, 
such as band pass filters and hig h Q re sonators, nonlin- 
earities are highly undesired [HJ, HH U2, 03 ■ 

Various nonlinear effects in superconductors and in 
NbN in particular have been reported and analyzed in the 
past by several research groups. Duffing like nonlinearity 
for example was observed in superconducting resonators 
employing different geometries and materials. It was ob- 
served in a high-T c superconducting (HTS) parallel plate 
resonator 01 , m a Nb microstrip resonator , in a Nb 
and NbN stripline resonators [l3|, in a YBCO coplanar- 
waveguide resonator 0] , in a YBCO thin film dielectric 
cavity 11 71 . and also in a suspended HTS thin film res- 
onator 18]. While other nonlinearities including notches, 
anomalies developing at the resonance lineshape and fre- 
quency hysteresis were reported in Refs. |2(J, |21] . 

However, in spite of the intensive study of nonlineari- 
ties in superconductors in the past decades, and the great 
progress achieved in this field, the determination of the 
underlying mechanisms responsible for microwave non- 
linear behavior of both low- and high-T c superconduc- 
tors is still a subject of debate This is partly be- 
cause of the variety of preparation and characterization 
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techniques employed, and the numerous fabrication pa- 
rameters involved. In addition nonlinear mechanisms in 
superconductors, which are usually divided into intrin- 
sic and extrinsic, are various and many times act con- 
currently. Thus identifying the dominant mechanism is 
generally indirect [23j . 

Among the nonlinear mechanisms investigated in su- 
perconductors one can name, Meissner effect |24j, pair- 
breaking, global and local heating effects 0, Il7| . rf and 
dc vortex penetration and motion |25| . defect points, 
damaged ed ges prjj. substrate material [27|], and weak 
links (WL) [28|]. Where WL is a collective term rep- 
resenting various material defects such as, weak super- 
conducting points switching to normal state under low 
current density, Josephson junctions forming inside the 
superconductor structure, grain-boundaries, voids, insu- 
lating oxides, insulating planes. These defects and im- 
purities generally affect the conduction properties of the 
superconductor and as a result cause extrinsic nonlinear 
effects. 

In this paper we report the observation of unique non- 
linear effects measured in the resonance lineshape of NbN 
superconducting microwave resonators. Among the ob- 
served effects, asymmetric resonances, multiple jumps in 
the resonance curve, hysteretic behavior in the vicinity of 
the jumps, frequency hysteresis loops changing direction, 
jump frequency shift as the input power is increased and 
nonlinear coupling. Some of these nonlinear effects were 
introduced by us in a previous publication |29|. Thus 
this paper will focus on presenting a new set of mea- 
surements applied to these nonlinear resonators, which 
provides a better understanding of the underlying physi- 
cal mechanism causing these effects. To this end, we have 
measured the nonlinear superconducting resonators using 
different operating conditions, such as bidirectional fre- 
quency sweeps, added white noise, fast frequency sweep 
using frequency modulation (FM), and dc magnetic field. 
In each case we observe a unique nonlinear dynamics of 
the resonance lineshape which is qualitatively different 
from the commonly reported Duffing oscillator nonlin- 
earity. These nonlinear effects are shown to originate 
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from WL forming at the boundaries of the NbN columnar 
structure. A theoretical model explaining the dynamical 
behavior of the resonance lineshape in terms of abrupt 
changes in the macroscopic parameters of the resonator 
is formulated. While these abrupt transitions in the char- 
acteristic parameters of the resonator are attributed to 
local heating of WL. Furthermore, simulations based on 
this model are shown to be in a very good qualitative 
agreement with the experimental results. 

The remainder of this paper is organized as follows, 
the fabrication process of the NbN superconducting res- 
onators is described briefly in Sec. II. The nonlinear re- 
sponse of these resonators measured using various operat- 
ing conditions are reviewed in Sec. III. Comparison with 
other nonlincarities reported in the literature is brought 
in Sec. IV. Possible underlying physical mechanisms re- 
sponsible for the observed effects are discussed in Sec. 
V. Whereas in Sec. VI, a theoretical model based on 
local heating of weak links is suggested, followed by sim- 
ulations qualitatively reproducing most of the nonlinear 
features observed in the experiments. Finally in Sec. VII, 
a short summary concludes this paper. 
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FIG. 1: (a) Schematic cross section of the stripline geometry 
used, consisting of five layers: two superconducting ground 
planes, two sapphire substrates and a NbN film in the middle 
deposited on one of the sapphires, (b) Top view of the three 
resonator layouts (Bl, B2, B3) which were deposited as the 
center layer. 



II. FABRICATION PROCESS 

The measurement results presented in this paper 
belong to three nonlinear NbN superconducting mi- 
crowave resonators. The resonators were fabricated us- 
ing stripline geometry, consisting of two superconducting 
ground planes, two sapphire substrates, and a center strip 
deposited in the middle (the deposition was done on one 
of the sapphire substrates). Fig. ^ shows a schematic 
diagram illustrating stripline geometry and a top view of 
the three resonator-layouts. We will refer to the three 
resonators in the text by the names Bl, B2 and B3 as 
defined in Fig. ^ The sapphire substrates dimensions 
used were 34mm X 30mm X 1mm, whereas the coupling 
gap between the resonators and their feedline was set 
to 0.4mm in Bl and B3 and 0.5mm in B2 resonator. 
The resonators were dc-magnetron sputtered in a mixed 
Ar/N2 atmosphere, near room temperature. The pat- 
terning was done using standard UV photolithography 
process, whereas the NbN etching was performed by Ar 
ion-milling. The sputtering parameters, fabrication de- 
tails, design considerations as well as physical properties 
of the NbN films can be found elsewhere 29] . The crit- 
ical temperature T c of Bl, B2 and B3 resonators were 
relatively low and equal to 10. 7K, 6.8K and 8.9K respec- 
tively. The thickness of the NbN resonators were 2200 A 
in Bl, 3000 A in B2, and 2000 A in B3 resonator. 



III. NONLINEAR RESONANCE RESPONSE 

In the following subsections, we present experimental 
results emphasizing the different aspects of the nonlin- 
ear response exhibited by Bl, B2 and B3 resonators. 



In subsection A, a resonance response measurement ob- 
tained while varying the input rf power is presented, 
showing an abrupt and low power onset of nonlinearity. 
In subsections B, C and D, representative experimental 
results measured while scanning the resonance response 
in the forward and backward directions are shown, ex- 
hibiting respectively, hysteresis loops changing direction, 
metastability and multiple jumps. Whereas in subsection 
E, the dependence of the resonance lineshape on applied 
dc magnetic field is examined, showing a change in the 
direction of the jump and a jump vanishing features. All 
measurements presented were performed at liquid helium 
temperature 4.2K, and the results were verified using two 
configurations, immersing in liquid helium and in vac- 
uum. 



A. Abrupt onset of nonlinearity 

In Fig. [2] we present a Sn parameter measurement 
of Bl first mode using a vector network analyzer. At 
low input powers, the resonance response lineshape is 
Lorentzian and symmetrical. As the input power is in- 
creased gradually in steps of 0.01 dBm, the resonance 
response changes dramatically and abruptly at about 
—28.04 dBm. It becomes extremely asymmetrical, and 
includes two abrupt jumps at both sides of the resonance 
lineshape. The magnitude of the jumps at some input 
powers, can be as high as 16 dB. As the input power is 
increased the resonance frequency is red shifted and the 
jump frequencies shift outwards away from the center fre- 
quency. Moreover at much higher powers not shown in 
the figure, the resonance curve becomes gradually shal- 
lower and broader in the frequency span. It is also worth- 
while noting that the intensive evolution of the resonance 
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Frequency [GHz] 

FIG. 2: Sn amplitude measurement of Bl resonator at its 
first mode. The measured resonance lineshapes are asymmet- 
rical and contain two abrupt jumps at the sides of the reso- 
nance. Moreover the jump frequencies shift outwards as the 
input power is increased. The measured resonance lineshapes 
were shifted vertically by a constant offset for clarity. 



lincshapc shown in Fig. Stakes place within only 1 dBm 
power range. 



B. Hysteretic behavior 

As the equation of motion for nonlinear systems be- 
comes multiple valued or lacks a steady state solution 
in some parameter domain, nonlinear systems tend to 
demonstrate hysteretic behavior with respect to that pa- 
rameter. 

Frequency hysteresis in the resonance lineshape of su- 
perconducting resonators exhibiting Duffing oscillator 
nonlinearity and other kinds of nonlinearities were ob- 
served by several groups 0, y(J]. Hysteretic behavior 
and losses in superconductors were discussed also in Refs. 
[25H3ll |. Moreover recent works, which examined the res- 
onance response of rf tank circuit coupled to a SQUID, 
have reported several interesting frequency hysteresis fea- 
tures [Him!!!. 

Likewise measuring the resonance response of our non- 
linear resonators yields a hysteretic behavior in the vicin- 
ity of the jumps. However, this hysteretic behavior 
is unique in many aspects. In Fig. we show a 
S11 measurement of Bl resonator at its first mode, mea- 
sured while sweeping the frequency in both directions. 
The input power range shown in this measurement corre- 
sponds to a higher power range than that of Fig. [21 The 
red line represents a forward frequency sweep, whereas 
the blue line represents a backward frequency sweep. 

At —20.6 dBm the resonance lineshape contains two 
jumps in each scan direction and two hysteresis loops. 
The left hysteresis loop circulates clockwise whereas the 
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FIG. 3: (Color online). Frequency sweep measurement of Bl 
resonator at its first mode performed in both frequency direc- 
tions. The plots exhibit hysteresis loops forming at the vicin- 
ity of the jumps, as well as hysteresis loops changing direc- 
tion as the input power is increased. The red line represents 
a forward sweep, whereas the blue line represents a backward 
sweep. The number of measurement points employed in each 
scan direction is 500 points. The resonance lineshapes were 
shifted vertically by a constant offset for clarity. In the inset, 
a "Zoom in" measurement of the right hysteresis loop of Bl 
first resonance is shown. The measurement, obtained using a 
spectrum analyzer, includes 100 data points and corresponds 
to —20.3 dBm input power. 



right loop circulates counter clockwise. However the com- 
mon property characterizing them is that the jumps oc- 
cur at higher frequencies in the forward scan compared 
to their counterparts in the backward scan. As the input 
power is increased to about —20.2 dBm the two opposed 
jumps at the left side meet and the left hysteresis loop 
vanishes. At about —19.4 dBm a similar effect happens to 
the right hysteresis loop, and it vanishes as well. Whereas 
at higher input powers (i.e. —19 dBm, —18.6 dBm) the 
two jumps occur earlier at each frequency sweep direc- 
tion, causing the hysteresis loops to appear circulating 
in the opposite direction compared to the —20.6 dBm 
resonance curve for instance. As we show in the next 
subsection, this picture of well defined hysteresis loops is 
strongly dependent on the applied frequency sweep rate 
and on the system noise. A possible explanation for this 
unique hysteretic behavior would be presented in Sec. 
VI. 



C. Metastable states 

Jumps in the resonance response of a nonlinear oscil- 
lator are usually described in terms of metastable/stable 
states and dynamic transition between basins of attrac- 
tion of the oscillator [3j|, thus in order to examine the 
stability of these observed resonance jumps, we carried 
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out several measurements. 

In one measurement, we have "zoomed in" around the 
right jump of the resonance at —20.3 dBm and exam- 
ined its frequency response in both directions. The mea- 
surement setup included a signal generator, the cooled 
resonator and a spectrum analyzer. The reflected signal 
power off the resonator was redirected by a circulator, 
and measured using a spectrum analyzer. The measure- 
ment result obtained using 100 sampling points in each 
direction, is exhibited in the inset of Fig. OH where the 
metastable nature of the jump region is clearly demon- 
strated. 

In another measurement configuration we have investi- 
gated this metastability further by monitoring the effect 
of applied broadband noise on the resonance jumps. We 
applied a constant white noise power to the resonator, 
several orders of magnitude lower than the main signal 
power, using the setup depicted in Fig. 21 The applied 
white noise level was —58 dBm/Hz (measured separately 
using spectrum analyzer), and was generated by amplify- 
ing the thermal noise of a room temperature 50 Q load us- 
ing an amplifying stage. The generated noise was added 
to the transmitted power of a network analyzer via a 
power combiner. The power reflections were redirected 
by a circulator, and were measured at the second port 
of the network analyzer. The effect of the —58 dBm/Hz 
white noise power on Bl first mode jumps, is shown in 
Fig. 03 (a), whereas in Fig. |S](b) we show for comparison 
the nearly noiseless case obtained after disconnecting the 
amplifier and the combiner stage. The two measurements 
were carried out in the same input power range (—23.9 
dBm through —20 dBm). 

By comparing between the two measurement results, 
one can make the following observations. The two fold 
jumps in Fig. 03(b) form a hysteresis loop at both sides of 
the resonance curve. By contrast in Fig. 03(a), as a result 
of the added noise, the hysteresis loops at the right side 
vanish, while the jumps at the left side, become frequent 
and bidirectional (indicated by the thick colored lines). 

At a given input drive, the transition rate T (/) 
between the oscillator basins of attraction, can 
be generally estimated by the expression T (f) — 
r exp (-Ea (/) /k B T eff ) PHI, where E A (/) is the quasi- 
activation energy of the oscillator, T e ff is proportional 
to the noise power, ks is Boltzmann's constant, / is the 
oscillator frequency, whereas Tq is related to Kramers 
low-dissipation form [3t1 | and it is given approximately 
by fo/Q, where fo is the natural resonance frequency, 
and Q is the quality factor of the oscillator. From our 
results we roughly estimate the order of magnitude of Ea 
to be 10 14 K for the jump on the left Note however 
that this quantity varies between different transitions and 
strongly depends on the operating point. 
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FIG. 4: Schematic diagram of the experimental setup used to 
measure the nonlinear first resonance of Bl using frequency 
sweep mode, while applying white noise. 
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FIG. 5: (Color online). Frequency sweep measurement of 
Bl resonator first mode performed in both directions while 
(a) applying white noise of —58 dBm/Hz (b) without apply- 
ing external noise. The red line represents a forward sweep, 
whereas the blue line represents a backward sweep. The mea- 
sured resonance curves were shifted vertically by a constant 
offset for clarity. 



D. Multiple jumps 

Another nonlinear feature, namely multiple jumps in 
the resonance lineshape, is observed when measuring the 
resonance response of B3, while sweeping the frequency 
in the forward and backward directions. In Fig. 03 wc 
show a representative measurement of the first resonance 
of B3 corresponding to 1.49 dBm input power, exhibiting 
three jumps in each sweep direction and four hysteresis 
loops. 



E. Dc magnetic field dependence 

Measuring B2 resonator second mode under dc mag- 
netic field yielded additional nonlinear features in the 
resonance response lineshape as shown in Figs. [7| 00 

In Fig. [7| we show the resonance lineshape of B2 sec- 
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FIG. 6: (Color online). Sn parameter amplitude measure- 
ment of the first resonance of B3 resonator, (shown in the 
inset), measured at input power of 1.49 dBm. The measure- 
ment was done using a network analyzer employing 4000 mea- 
surement points in each direction. The red line represents 
a forward frequency scan whereas the blue line represents a 
backward scan. The plot shows clearly three jumps within the 
resonance lineshape in each direction, as indicated by small 
circles. 



ond mode measured while applying a perpendicular dc 
magnetic field of 90mT. As the input rf power is in- 
creased gradually, the resonance lineshape undergoes dif- 
ferent phases. While at low and high powers the curves 
are Lorentzians and symmetrical, in the intermediate 
range, the resonance curves include a jump at the left 
side, which, as the input power increases, flips from the 
upward to the downward direction. 

Whereas in Fig. [3J where we have set a constant input 
power of —5 dBm and increased the applied magnetic 
field by small steps, the left side jump vanishes as the 
magnetic field exceeds some relatively low threshold of 
~ 11.8mT. These effects will be further discussed in 
Sees. V and VI. 



IV. COMPARISON WITH OTHER 
NONLINEARITIES 
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FIG. 7: B2 nonlinear resonance response measured under a 
constant magnetic field of 90mT, while increasing the input 
power. The resonance which starts as a Lorentzian at low 
powers, develops into a resonance curve having an upward 
jump, a curve with no jump, a curve having a downward 
jump and finally a Lorentzian curve again as the power is 
increased. The measured curves were shifted vertically by a 
constant offset for clarity. 
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FIG. 8: Increasing the magnetic field gradually from zero 
causes the jump in the B2 resonance lineshape to disappear 
at relatively low value of 11.8mT, while applying a constant 
input power level of —5 dBm. This jump vanishing indicates 
that the jump mechanism is sensitive to magnetic field. The 
measured curves were shifted vertically by a constant offset 
for clarity. 



The most commonly reported nonlinearity in super- 
conductors is the Duffing oscillator nonlinearity. How- 
ever this nonlinearity is qualitatively different from the 
nonlinearity we observe in our NbN samples and which 
is reported in this paper. In Fig. we show, for the 
sake of visual comparison, a resonance response measured 
at 4.2K, exhibiting Duffing oscillator nonlinearity of the 
kind generally reported in the literature p^ . ll4llTq . This 
nonlinearity which can be explained in terms of resistance 
change AR and kinetic inductance change ALk 0, H(j 
was measured at the first resonance frequency of a 2200 



A thickness Nb resonator employing B2 layout geometry 
(T c = 8.9K). The differences between the two nonlinear 
dynamics shown in Figs. |3 and ^\ are obvious. In Fig. 
[5] the nonlinearity is gradual, while in Fig. [2 the power 
onset of nonlinearity is abrupt and sudden. In Fig. 
the resonance response in the nonlinear regime, contains 
an infinite slope at the left side, whereas in Fig. the 
curves contain two jumps at both sides of the resonance 
response. In Fig. [§] changes in the resonance curve are 
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FIG. 9: Duffing oscillator nonlinearity exhibited by a Nb 
stripline resonator employing B2 layout at its first mode. The 
different Sn amplitude plots correspond to different input 
powers, ranging from —15 dBm to 15 dBm in steps of 1 dBm. 
As the input power is increased the resonance becomes asym- 
metrical and infinite slope builds up at the left side of the 
resonance curve. The plots were offset in the vertical direc- 
tion for clarity. 



measured on power scale of 1 dBm, while changes in Fig. 
|2]are measured on 0.01 dBm scale. Whereas the onset of 
nonlinearity in Fig. [5]is of the order of 10 dBm, the onset 
of nonlinearity in Fig. [21 is about 4 orders of magnitudes 
lower ~ —28 dBm. Furthermore, by taking into account 
the differences in the hysteretic behavior of the two non- 
linearities and the multiple jumps feature shown in Fig. 
[SI the special characteristics of the reported nonlinearity 
are further established. 

Abrupt jumps in the resonance lineshape similar in 
some aspects to the jumps reported hereim were observed 
in two-port high-T c YBCO resonators [HIHEI]. Portcs 
et al. |l!j have also reported some frequency hysteretic 
behavior in the vicinity of the jumps. However, one sig- 
nificant difference between the two nonlinearities is the 
onset power of nonlinearity reported in these references, 
which is on the order of 20 dBm [2(], , that is about 5 
orders of magnitude higher than the onset power of non- 
linearity of Bl first mode (~ —28 dBm). All three works 
O I20I l2l| have attributed the nonlinear abrupt jumps 
to local heating of distributed WL in the resonator film. 



V. POSSIBLE NONLINEAR MECHANISMS 



producible results with a good accuracy in the magnetic 
field magnitude, the microwave input power and in the 
jump frequency (less than 200 kHz offset). Moreover the 
low magnetic field threshold ~ 11.8 mT above which the 
B2 resonance jump vanished, is about 3.5 times lower 
than He, (flux penetration) of NbN reported for exam- 
ple in jl3| . 

In the following subsection A, we provide a direct evi- 
dence of WL. Whereas in subsection B, we exclude global 
heating mechanism as a possible source of the effects. 



Columnar structure 



It is well known from numerous research works done in 
the past [4ll l42l | , that NbN films can grow in a granular 
columnar structure under certain deposition conditions. 
Such columnar structure may even promote the growth 
of random WL at the grain boundaries of the NbN films. 
To study the NbN granular structure we have sputtered 
about 2200A NbN film on a thin small rectangular sap- 
phire substrate of 0.2 mm thickness. The sputtering con- 
ditions applied were similar to those used in the fabrica- 
tion of B2 resonator. Following the sputtering process, 
the thin sapphire was cleaved, and a scanning electron 
microscope (SEM) micrograph was taken at the cleavage 
plane. The SEM micrograph in Fig. ^3 which clearly 
shows the columnar structure of the deposited NbN film 
and its grain boundaries, further supports our weak link 
hypothesis. The typical diameter of each NbN column is 
about 20 nm. 
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The relatively very low onset power of nonlinearity ob- 
served in these resonators as well as its strong sensitivity 
to rf power, highly imply an extrinsic origin of these ef- 
fects, and as such, hot spots in WL is a leading candidate 
for explaining the nonlinearity. 

Vortex penetration in the bulk or WL is less likely, 
mainly because heating the sample above T c between se- 
quential magnetic field measurements, have yielded re- 



FIG. 10: A SEM micrograph displaying a 2200A NbN film de- 
posited on a thin sapphire substrate using similar sputtering 
conditions as B2 resonator. The micrograph exhibits clearly 
the columnar structure of the NbN film and its grain bound- 



B. Frequency sweep time analysis 

Resistive losses and heating effects are typically char- 
acterized by relatively long time scales [l4(- In attempt 
to consider whether such effects are responsible for the 
observed nonlinearities in general and for the jumps in 
particular, we have run frequency sweep time analysis 
using the experimental setup depicted in Fig. ^2 We 
have controlled the frequency sweep cycle of a signal gen- 
erator via FM modulation. The FM modulation was ob- 
tained by feeding the signal generator with a saw tooth 
waveform having 1// sweep time cycle. The reflected 
power from the resonator was redirected using a circula- 
tor and measured by a power diode and oscilloscope. The 
left and right hand jumps of B2 ~ 4.39GHz resonance 
were measured using this setup, while applying increas- 
ing FM modulation frequencies up to 200kHz. In Fig. 
1121 we present a measurement result obtained at 50kHz 
FM modulation, or alternatively T sweep of 20/iS. The FM 
modulation applied was ±20MHz around 4. 4022GHz cen- 
ter frequency. The measured resonance response appears 
inverted in the figure due to the negative output polarity 
of the power diode. The fact that both jumps continue 
to occur within the resonance lineshape (see Fig. I12fl . 
in spite of the short duty cycles that are of the order of 
~ /is, indicates that heating processes which have typi- 
cal time scale on the order of s to ms 14j are unlikely to 
cause these effects. 
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FIG. 11: Frequency sweep time analysis setup. The frequency 
sweep time of the microwave signal generator was FM mod- 
ulated by a saw tooth waveform with frequency /. The re- 
flected power from the resonator was measured by a power 
diode and oscilloscope. 




FIG. 12: Frequency sweep time measurement. The figure 
displays the resonance measured by Agilent oscilloscope while 
applying a saw tooth FM modulation of frequency 50kHz ( 
Tsweep =20/is) to the signal generator. The left and right 
jumps of the resonance are still apparent in spite of the fast 
rate frequency sweep. Thus indicating that the jumps do not 
originate from any global heating mechanism. 



where the parameters C ~ 2.7 • 10 _3 Jcm _3 K _:l (NbN) 
H3, d = 3000A (B2 thickness), and a ~ 1.5Wcm- 2 K-i 
at 4.2K (sapphire substrate) [40j . have been used. Sim- 
ilar calculation based on values given in Ref. [3^| yields 
t ~ 2.1 • 10~ 9 s. These time scales are of course 2-3 orders 
of magnitude lower than the time scales examined by the 
FM modulation setup, and thus local heating of WL is 
not ruled out. 



VI. LOCAL HEATING OF WL MODEL 

In this section we consider a hypothesis according to 
which local heating of WL is responsible for the observed 
effects. We show that this hypothesis can account for the 
main nonlinear features observed, and that simulations 
based on such a theoretical model, exhibit a very good 
qualitative agreement with experimental results. 



A. Theoretical modeling 



However the above measurement result does not ex- 
clude local heating of WL [H El E3|- Assuming that 
the substrate is isothermal and that the hot spot is dissi- 
pated mainly down into the substrate rather than along 
the film [23, one can evaluate the characteristic relax- 
ation time of the hot spot using the equation r = Cd/a, 
where C is the heat capacity of the superconducting film 
(per unit volume), d is the film thickness, and a is the 
thermal surface conductance between the film and the 
substrate [sH- Substituting for our B2 NbN resonator 
yields a characteristic relaxation time of r ~ 5.4 • 10 _8 s, 



Consider a resonator driven by a weakly coupled feed- 
line carrying an incident coherent tone b tn e~ iupt , where 
b ln is a constant complex amplitude and lo p is the drive 
angular frequency. The mode amplitude inside the res- 
onator A can be written as A = Be -lWpt , where B (t) is 
a complex amplitude, which is assumed to vary slowly 
on the time scale of 1/lu p . In this approximation, the 
equation of motion of B reads 



— = li(Lo p -o Jo )-j}B~i^- 1 b m + c m , (1) 
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where loq is the angular resonance frequency, 7 = 71 + 
72 , 71 is the coupling constant between the resonator and 
the feedline, and 72 is the damping rate of the mode. The 
term c m represents input noise with vanishing average 



temperature of other parts of the resonator is assumed 
be equal to that of the coolant T . The power Q heating 
up the hot spot is given by Q = aQt where ^ a ^ 1. 
The heat balance equation reads 



(O = 0, 

and correlation function given by 



{c m {t)c m *(t')) =Gcj S{t-t') 



(2) 



(3) 



In thermal equilibrium and for the case of high temper- 
ature kgT >> hu>o, where kg is Boltzmann's constant, 
one has 



G 



27 k B T 

UJQ flUJQ 



(4) 



In terms of the dimensionless time r = uj t, Eq. 
reads 



(10) 



where C is the thermal heat capacity, W = H (T — To) 
is the power of heat transfer to the coolant, and H the 
heat transfer coefficient. Defining the dimensionless tem- 
perature [H ] 



e 



T-Tn 



T c - To ' 

where T c is the critical temperature, one has 



(11) 



d® 

~dw 



-5(© -0c 



where 



(12) 



f= li ^- UJ ° ) -" (B-B oo)+ C -, (5) 
dr ujq ujq 



where 



i (uj p - w ) - 7 



(6) 



Small noise gives rise to fluctuations around the steady 
state solution B^. A straightforward calculation yields 



B-B 



Gujq 
2 7 ' 



(7) 



be 



The output signal a out reflected off the resonator can uc 
written as a ouf = & 0, "t e -«*>pt. The input-output relation 
relating the output signal to the input signal is given by 



c 



p = 



^0.9 



C(T C -T ) 



H 



(13) 



(14) 



(15) 



While in Duffing oscillator equation discussed in Ref. 
[2| , the nonlinearity can be described in terms of a grad- 
ually varying resonance frequency dependent on the am- 
plitude of the oscillations inside the cavity, in the current 
case, the resonance frequency ujo, the damping rates 71, 
72 and a factor are considered to have a step function 
dependence on T, the temperature of the WL 




(8) 



Whereas the total power dissipated in the resonator Q t 
can be expressed as yj 



UJ 



lo 0s e < 1 

W 0n > 1 



71, © < 1 

7m © > 1 ' 



(16) 



(17) 



(9) 



where E = \B\ . 

Furthermore, consider the case where the nonlinear- 
ity is originated by a local hot spot in the stripline res- 
onator. If the hot spot is assumed to be sufficiently small, 
its temperature T can be considered homogeneous. The 



72 



72s © < 1 
72n © > 1 ' 



a s 6 < 1 
a n 6 > 1 



(18) 



(19) 
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In general, while disregarding noise, the coupled dif- 
ferential equations (JSJ and (|T2|l may have up to two dif- 
ferent steady state solutions. A superconducting steady 
state of the WL exists when Ooo < 1, or alternatively 
when E < E s , where E s — gC (T c — To) j2a s ^i s h. Sim- 
ilarly, a normal steady state of the WL exists when 
0oo > 1, or alternatively when E > E n , where E n = 
gC (T c - T ) /2a„72„fi. 

In general, the reflection coefficient Su in steady state 
is given by [3| 



Su 



^out 

b in 



72 - 7i - i {Up - up) 
72 + 7i - i (w P - w ) ' 



(20) 



B. Simulation results 



Simulating the resonator system using this local heat- 
ing WL model, yields results which qualitatively agree 
with most of the nonlinear effects previously presented. 
In subsection (1) we simulate the main effects of Sec. Ill 
(A, B, C, D), whereas in subsection (2) we simulate and 
provide a possible explanation to the nonlinear features 
of Sec. Ill (E). 



1. Abrupt jumps and hysteretic behavior 

In Fig. 1131 we show a resonance response simulation 
result based on the hot spot model, which simulates the 
abrupt jump exhibited in Fig. [21 The solid and dotted 
lines represent valid steady state solutions of the system 
and invalid steady state solutions respectively. Whereas 
the blue and red colors represent superconducting WL 
solutions and normal WL solutions respectively. In plot 
(a) the superconducting WL solution is valid in the nor- 
malized frequency span, and therefore the system follows 
this lineshape without jumps. As we increase the am- 
plitude drive bi n we obtain a result seen in plot (b). As 
the frequency is swept , jumps in the resonance response 
are expected to take place as the solution followed by 
the system (according to the initial conditions) becomes 
invalid. Thus in the forward sweep direction (as the fre- 
quency sweep of Fig. we get two jumps indicated by 
black arrows on the figure. Similar to Fig. |3 the mag- 
nitudes of the jumps in the plot are unequal (the left 
jump is higher). This difference in the magnitude of the 
jumps is generally dependent on the relative position be- 
tween the two resonance frequencies (Eq. 116(1 . while in 
measurement due to the metastability of the system in 
the hysteretic regime, it depends also on the frequency 
sweep rate. The simulation parameters used in the dif- 
ferent cases are indicated in the figure captions. 

The behavior of the frequency hysteresis loops exhib- 
ited in Fig. |2| is simulated in Fig. The different 
plots exhibited in Fig. ^] correspond to the different 
cases shown in Fig. [3J In plot (a) the jumps in the for- 
ward direction (indicated by the arrows in that direction) 



b. =83900 



_ -30 

g 0.99 0.992 0.994 0.996 0.998 
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1.002 1.004 1.006 1.008 1.01 
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FIG. 13: (Color online). Simulated resonance response ob- 
tained by the hot spot mechanism modeling. The differ- 
ent plots simulate the nonlinear behavior shown in Fig. [5] 
Panels (a) and (b) correspond to an increasing drive am- 
plitude b m . The solid lines represent valid steady state so- 
lutions whereas the dotted lines represent invalid solutions. 
The red lines represent the normal WL solutions, whereas the 
blue lines represent the superconducting WL solutions. The 
black arrows show the direction of the jumps in the different 
cases. The parameters that were used in the simulation are: 
tW^On = 1.0008, 7i„/w n = 2.5- 10" 3 , -y ls /Luo n = 1.5- 10" 3 , 
T2„K = 2.75 • 10~ 3 , 7 2s /a>on = 5 ■ 10" 3 , a n = 0.8, a s = 1, 
g = 0.5, p= 10" 10 . 



occur at higher frequencies than the jumps in the back- 
ward direction. Whereas in plot (b) corresponding to a 
higher amplitude drive bi n we show a case in which the 
left side hysteresis loop vanish as the two opposed jump 
frequencies coincide. If we increase bt n further, then at 
some amplitude drive level as shown in plot (c), we get a 
similar case of hysteresis loop vanishing at the right side 
of the resonance response. Whereas at the left side we 
get a frequency region where both the superconducting 
and the normal WL solutions are invalid. In this instable 
region, transitions between the invalid solutions are ex- 
pected, depending on the number of the sampling points, 
the sweep time, and the internal noise. However due to 
this instability, the system is highly expected to jump 
"early" in each frequency direction, as it enters this re- 
gion (at lower frequencies in the forward direction, and 
at higher frequencies in the backward direction). Thus 
leading to the observed change in the direction of the 
hysteresis loop. By increasing 6j„ further, one obtains a 
case in which both hysteresis loops are circulating in the 
opposite direction compared to plot (a). 

Furthermore the intermediate jump indicating insta- 
bility, which appear at the left jump region of the last 
resonance curve in Fig. [3] (corresponding to —18.6 dBm), 
can be explained by this model as well. By solving the 
coupled Eqs. (J5J, (|12f> in the time domain for a single 
normalized frequency lo p = 0.9952 (arbitrarily chosen in 
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FIG. 14: (Color online). Simulated resonance response ob- 
tained by the hot spot mechanism modeling. Panels (a), (b), 
(c) and (d) correspond to an increasing drive amplitude b m . 
The different plots simulate the nonlinear behavior shown 
in Fig. [3] The lines and symbols are the same as in Fig. 
1131 The parameters that were used in the simulation are: 
u 0s /cj 0n = 0.99989, yi n /uon = 2.5-1CT 3 , y ls /u> 0n = 1.5-HT 3 , 
72n/w 0n = 2.75 • 10" 3 , 72s/won = 5 ■ 10" 3 , a„ = 0.8, a s = 1, 

-10 



3 = 0.5, p = 10" 
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FIG. 15: The dynamical solution of the coupled Eqs. 1121 
at a normalized frequency of u v — 0.9952. The simulation 
parameters used are the same as those of plot (d) in Fig. 1141 



the left side hysteresis region) and using the simulation 
parameters of plot (d) in Fig. El one obtains the oscilla- 
tion pattern of the dimcnsionless parameter shown in 
Fig. El as a function of the dimcnsionless time r. The 
oscillations indicating instability are between the su- 
perconducting and the normal values, corresponding to 
6 < 1 and 6 > 1 respectively. 

As to the multiple jumps feature exhibited in Fig. |U 
a straightforward generalization of the model may be 
needed in order to account for this effect. Such gener- 
alization would include WL having a variation in their 
sizes and critical current along the stripline. Thus caus- 



ing them to switch to normal state at different current 
drives (at different frequencies) and as a result induce 
more than two jumps in the resonance lincshape. 



2. Magnetic field dependence 

In this subsection we show how the model of local heat- 
ing of WL can also account for the nonlinear dynamics of 
the resonance lineshape observed under applied magnetic 
field. 

To this end, we show in Fig. 1161 a simulation result 
based on the WL local heating model, which qualitatively 
regenerates the nonlinear behavior of the resonance line- 
shape of B2 under a constant magnetic field (presented 
in Fig. UJ. At low drive amplitude bi n , only the super- 
conducting WL steady state solution exists, and thus no 
jump occurs as one sweeps the frequency (plot (a)). In- 
creasing the drive amplitude bi n (plot (b)) causes the su- 
perconducting WL solution to become invalid in the cen- 
ter frequency region, thus the resonance response jumps 
upward (as the system reaches the invalid region) and 
stabilizes on the normal WL steady state solution, as in- 
dicated by arrows on the plot. By increasing the drive 
amplitude further (plot (c)) there exists an intersection 
point where a smooth transition without a jump is ex- 
pected to occur between the valid superconducting WL 
solution and the valid normal one. Whereas in plot (d) 
where we have increased bi„ further, a downward jump 
in the resonance response occurs as the valid normal WL 
solution lies below the invalid superconducting WL solu- 
tion. Finally in plot (e) corresponding to a much higher 
drive, only the normal WL steady state solution exists 
within the frequency span and therefore there are no 
jumps in the resultant curve. 

Another measurement which can be explained using 
the WL model, is the measurement shown in Fig. [SJ 
where the left side jump vanishes as the magnetic field 
increases above some low magnetic field threshold. This 
result can be explained in the following manner. Increas- 
ing the applied dc magnetic field would inevitably in- 
crease the screening supercurrent flowing in the film and 
the local heating of WL. As the local heating exceeds 
some threshold, the superconducting WL solution would 
become invalid (for the same frequency span), and con- 
sequently the system would only follow the normal WL 
solution without apparent jumps. 



VII. SUMMARY 

In attempt to investigate and manifest nonlinear effects 
in superconducting microwave resonators, several super- 
conducting NbN resonators employing different layouts, 
but similar sputtering conditions have been designed and 
fabricated. The resonance lineshapes of these NbN res- 
onators having very low onset of nonlinearity, several 
orders of magnitude lower than other reported nonlin- 
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FIG. 16: (Color online). Simulated resonance response ob- 
tained by the hot spot mechanism modeling. The different 
plots simulate the nonlinear behavior shown in Fig. [7| Pan- 
els (a), (b), (c), (d) and (e) correspond to an increasing drive 
amplitude b m . The lines and symbols are the same as in Fig. 
1181 The simulation parameters used are: uio n /uio a = 1.0002, 
y ln /u>0s = 0.0029, yi s /u 0s = 0.0019, 7 2 „/^ 0s = 0.0019, 
72,/W = 0.0015, a„ = 1, a s = 0.8, g = 0.5, p = KT 10 . 



earities 0, OH C3 > exhibit some extraordinary nonlin- 
ear dynamics. Among the nonlinearities observed, while 
applying different measurement configurations, abrupt 
metastable jumps in the resonance lineshape, hystere- 



sis loops changing direction, multiple jumps, vanishing 
jumps and jumps changing direction. These effects are 
hypothesized to originate from weak links located at the 
boundaries of the columnar structure of the NbN films. 
This hypothesis is fully consistent with SEM micrographs 
of these films, and generally agrees with the extrinsic like 
behavior of these resonators. To account for the various 
nonlinearities observed, a theoretical model assuming lo- 
cal heating of weak links is suggested. Furthermore, sim- 
ulation results employing this model are shown to be in 
a very good qualitative agreement with measurements. 

Such strong sensitive nonlinear effects reported herein 
may be utilized in the future in a variety of applications, 
ranging from qubit coupling in quantum computation, 
to signal amplification [43 and to the demonstration of 
some important quantum effects in the microwave regime 

urn. 
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